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I. INTRODUCTION 



Peter H0yer* 

BRICS^ , Department of Computer Science, University of Aarhus, 
Ny Munkegade, Bldg. 540, DK-8000 Aarhus C, Denmark. 

We consider the use of arbitrary phases in quantum amplitude amplification which is a generali- 
zation of quantum searching. We prove that the phase condition in amplitude amplification is given 
by tan(y/2) = tan(0/2)(l — 2a), where cf> and cp are the phases used and where a is the success 
probability of the given algorithm. Thus the choice of phases depends nontrivially and nonlinear ly 
on the success probability. Utilizing this condition, we give methods for constructing quantum 
algorithms that succeed with certainty and for implementing arbitrary rotations. We also conclude 
that phase errors of order up to can be tolerated in amplitude amplification. 
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Most quantum algorithms developed so far are based on two techniques: Quantum Fourier Transforms and Ampli- 
7-H ' tude Amplification. The latter technique is a generalization of Grover's quantum algorithm for searching an unordered 
| database Q|, and it allows a quadratic speedup over any classical algorithm for many computational problems. Since 
i amplitude amplification is fundamental for quantum algorithms it has received a great deal of attention. This includes 
a study of its robustness to errors and modifications. In particular, the effects of using arbitrary phases in amplitude 
amplification have been studied in a sequence of papers by Long et al. ■ 

In this paper, we also consider the question of when arbitrary phases can be utilized in amplitude amplification. Our 
results complement the results of Long et al. whom are primarily interested in the question of how large phase errors we 
can tolerate and still obtain a quantum algorithm for searching that succeeds with high probability. We are primarily 
interested in the question of what restrictions we need to put on the two angles used in amplitude amplification and 
still obtain quantum algorithms that succeed with certainty. 

To illustrate our results, consider a database of size N with a unique marked element. In each iteration of Grover's 
algorithm, we rotate the phase of some states by angles <f> and <p, respectively. The main result in is that if the 
two angles differ by at most ^7= for some appropriate constant c, that is, if \ip — <f>\ < -^7=, then we can still find the 

marked element with high probability using only Q(^/N) iterations. Our main result is that if the two angles satisfy 
the equation tan(<p/2) = tan(</>/2)(I — j|) then we can find the marked element with certainty using only Q(y/N) 
iterations. Together with the result of Long et al., this provides a description of the use of arbitrary phases for 
bounded-error and exact quantum algorithms. It is possible to rederive the main result of Long et al. from our results 
by considering the case cj) = if an approximation to the perfect case tan(>/2) = tan(>/2)(I - -|). (See Sees. [Hi] 
below for rigorous statements.) 

We thus prove that there is phase condition in amplitude amplification, and that this condition is not that the 
phases are equal, but that they satisfy the trigonometric equation mentioned above. We believe that our approach is 
intuitive and that it yields short and straightforward proofs. 



II. AMPLITUDE AMPLIFICATION 



Amplitude amplification is a generalization of Grover's quantum searching algorithm that allows a speed up of 
many classical algorithms. The heart of amplitude amplification is an operator Q defined similarly to the operator 
used in Grover's algorithm fjl]. We refer the reader to IJ and the references therein for a throughly introduction to 
amplitude amplification. Here we give only a concise description of all the objects we require. 



* email: hoyer@brics . dk. 

' Basic Research in Computer Science, Centre of the Danish National Research Foundation. 



f 



Let TC be Hilbert space of dimension N and let {|0), ... , \N — 1)} be an orthonormal basis for 7i. Let A be a unitary 
operator on H. We may think of A as a quantum algorithm that uses no measurements. Let x '■ {0, • • • , N— 1} — > {0, 1} 
be a Boolean function. We say that a basis state |a;) is good if x( x ) — L an d otherwise we say that |a;) is bad. Given 
two angles < <f), ip < 2ir, define 



Q = Q(A, X , 4>, 9) = -ASo^A- 1 S x (<p). 
Here, the operator S x (ip) conditionally changes the phase of the amplitudes of the good states, 



(1) 



\X) 

\x) 



if xO) = 1 

if x(^) = o. 



(2) 



Similarly, the operator So(^) multiplies the amplitude by a factor of e 1 ^ if and only if the state is the zero state |0). 
Here and elsewhere we use i to denote the principal square root of —1. 

Let I'S) = -4|0) denote the superposition obtained by applying algorithm A on the initial state |0). Let = 



|*) where P goo d = 



x:xO) = l 



(x\ denotes the projection onto the subspace spanned by the good basis states, 
and similarly let |\I>o) = Pbadl*) where Pbad = J2 x - x ( x ) = o \ x )( x \- Let a = (4 f i|\I r i) denote the probability that a 
measurement of 1^) = _4|0) yields a good state, and let b = (^ol^o) denote the probability that a measurement of 
|*) yields a bad state. We then have that |*) = + |* ) and 1 = a + b. Finally, let angle 9 be so that < 9 < tt/2 
and a = sin 2 (0). 

As an example, for N = 2™, we obtain Grover's searching algorithm jl| by setting A to be the Walsh-Hadamard 
transform on n qubits, letting x( x ) be 1 if and only if the database holds a 1 at position x, and picking phases 
<f> = (p = 7r. If the database contains a 1 at t different positions then a — t/N . 

The operator Q implements a unitary operation on the subspace spanned by I'J'i) and l^o)- This subspace has 
dimension 2 if < a < 1. With respect to the ordered orthonormal basis ("^l^i): T^l^o)) > we can represent 
Q = Q(A, x, 4>, <p) by the 2x2 unitary matrix 



M = 



-{(l-e^)a + e^} (1 - e^^/aVT^ae 1 ^ 
(1 - e^vWT - ^ {(1 - e^)a - l}e ,v 



If (f> = ip = 7T, then this simplifies to Q| 



M 



cos(26>) -sin(26») 
sin(26») cos(26») 



(3) 



(4) 



That is, if we pick <f> — Lp — 7r, then each application of Q implements a rotation by angle 29. A natural question then 
is, what happens if at least one of the two angles <j> and ip is not equal to it? 



III. ARBITRARY ROTATIONS 

Consider the matrix M defined by Eq. (||). Our primary objective is to ensure that the diagonal elements of M are 
equal. 

Theorem 1 (Optimal angles) Suppose <p ^= it. Then the two diagonal elements of M are equal if and only if 

tan(^/2) = tan(>/2) (1 - 2a) (5) 

where matrix M is defined by Eq. (Qj. 

Equation (||) expresses the phase condition we want to impose on <fi and ip. Theorem 1 can be proven straight- 
forwardly using standard identities of the trigonometric functions. Figure |l| illustrates the proof for the case that 
< <j) < ir and < a < |. The other cases are proven similarly. 
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FIG. 1. Angle 95 as a function of 4> and a - Let £ be the length of the dotted line. Then tan(y>/2) = ( i — a) 4 and tan(0/2) ~ \ j, 
and hence tan(<^>/2) = tan(</!>/2)(l — 2a). 



We now show how to implement arbitrary rotations (up to certain phase factors) by appropriate choices of <fi and ip. 
Let < 1) < 2ir be any angle for which | sm(?9)| < sm(29). First we pick angle cf> so that the absolute value of the 
lower left entry of M equals | sin(?9)|. Then we pick angle p so that Eq. (JsJ) holds. This ensures that the two diagonal 
elements of M are equal. With these choices of <p and (p, matrix M can thus be written in the form 



H 



cos($) — sin(z?) 
sin($) cos(z9) 



1 



(G) 



for some angles < u, v < 2tt. We denote this matrix by H. Here and elsewhere, missing matrix entries are assumed 
equal to 0. 

To summarize, we have just shown that for all a (0 < a < 1) and for all angles t9 (0 < i9 < 2ir) so that | sin(z?)| < 
sin(26*) there exist angles < <fi, <p < 2tt so that M = H for some angles < u, v < 2ir, where M is given by Eq. (|^) 
and H by Eq. (^). Matrix H represents a rotation by angle 1? which is conjugated by a conditional phase change by 
angle u, up to a global phase factor of e m . 

In many applications, applying H will be equally good to applying a rotation by angle For instance, we can use H 
to implement arbitrary rotations as follows. Suppose a is known. Let x be any angle (0 < x < 2n). We implement 
a rotation by angle x as follows. First, we check if a: is a multiple of 29. If so, we simply just apply Q(-4, x, 7r, tt) a 
total number of x/(29) times and stop 0. Otherwise, we compute the smallest integer m larger than x/{29) and we 
set -d — x/m. Then we find angles </> and <p so that M = H, and we compute the angles u and v. With these choices 
of angles, we can factorize the rotation by angle x as 



cos(a:) 
sin(x) 



- sin(a;) 
cos(x) 



1 



M r 



1 



(7) 



Thus, to implement a rotation by angle x, we first apply a conditional phase-change of the bad states by angle u. Then 
we apply Q(.A, x,</>, p) a total number of m times, apply a conditional phase-change of the bad states by angle —u, 
and finally apply a global phase-change by angle -mn. 



IV. OBTAINING SUCCESS PROBABILITY 1 

We now show that if we pick nontrivial angles <p and tp so that Eq. (|^) holds, then we can find a good solution with 
certainty. That is, we can use any set of angles (</>, p) ^ (0, 0) for which Eq. (||) holds. 

Let angle < 1) < ir/2 be defined so that sin($) = | sin(0/2) sm(26*)|. Since </> and ip satisfy Eq. (||), we can write 
Q(-4, X, 4>, p) in the form 



"1 




cos($) 


- sin(tf)" 




"1 


e m 




sin(i?) 


cos($) 







(8) 



for some angles < u, v < 2n. That is, operator Q(-4, x, cf>, p) implements a rotation by angle i9, up to phase factors. 

Our idea for finding a good solution with certainty is as follows: Let m = [ 7r ^~ 9 "| and 9- ln \t = § — m&. Then 
— 9 < #j n it < 9. We first set up a superposition representing the initial angle 0jnjt- This is possible since |#i n it| < 0. 
Then we apply operator Q(^4, x,</>, <p) a total number of m times and finally we measure. This produces a good 
solution with certainty. 

A realization of this idea is as follows: Let l^init) denote the state sin(#i n i t )-^| , I , i) + cos(6>i n it) v / 1 1 _ a l^o)- First we 
apply A on the first register of the state |0)|0)|0), producing the state |^)|0)|0). Then we unitarily map |\I>)|0)|0) to 
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the state a|4'init)|0)|0)+/3|0)|^'i)|l) for some complex numbers a, (3 S C. Then we apply operator Q(A, x, (p, (p) atotal 
number of m times on the first register, producing the state a'| 1 E r i)|0)|0) + /3 / |£')|\f , i)|l) for some complex numbers 
a' , j3' € C and some state \E). Finally, we swap the contents of the first two registers conditionally to that the third 
register contains a 1, producing the final tensor product state |4 , i)|Z?') 2 3 , for some state \E') 2 3 that represents the 
joint state of registers 2 and 3. 

We may summarize this section by saying that for fixed known rotational angle d, we can modify the angle of the 
initial state so that we succeed with certainty. 



V. THE PHASE CONDITION 



Our condition that tan(</?/2) = tan(</>/2)(l — 2a) implies unfortunately that 4> and </S> depend nontrivially on a. Put 
formally, for all angles < <f> < 2tt so that <fi ^ tt, the following holds: For all angles < p < 2tt with p ^ tt, there 
exists a unique a £ S so that Eq. ^ holds, and for all < a < 1, there exists a unique angle < p < 2tt so that 
Eq. (|^) holds. If the success probability a is known in advance, then for any angle (f> we may want to pick, we can 
easily compute the angle <p to use. However, if a is not known, then this is not possible unless (f> <E {0,ir}. Thus for 
angles 4> ^ {0, tt}, we need to know a to compute <p so that Eq. (JsJ) holds. 

If we do not know a in advance, then a subsidiary strategy could be to utilize a set of angles (0, p) so that 
Eq. (§) almost holds. If a is small then we could for example approximate p by (f>. In the papers []2H5|, Long et al. 
consider the question of when arbitrary phases can be utilized successfully in quantum searching. Their conclusion 
is that the angles have to equal (<fr — <p) "to construct an efficient quantum search algorithm" ||. Our condition 
that tan(</?/2) = tan(</>/2)(l — 2a) is obviously different from their condition (that </> = ip) whenever < a < 1 and 
<j) £ {0, 7r}. The explanation for these different results is that Long et al. consider when the quantum search algorithm 
succeeds with high probability, whereas we, in the previous section, consider when the quantum search algorithm 
succeeds with certainty. In particular, all our calculations are exact. A main proof technical idea used by Long et al. 
for example in pf is approximations of the type Ki « e K2 for 2 x 2 matrices Ki and 

Long et al. |pH prove their result via an SO(3) rotational interpretation of operator Q(^4, X, </>, <p). We now reprove 
the theorem of Long et al. that if we use phases <f> = <p, then we can find a good solution with high probability. The 
main idea in our alternative proof is to consider the case <fi — p an approximation to the perfect case in which Eq. (§) 
holds, and then lower bound how well this approximation works. We do that by upper bounding the norm of the 
difference of two operators. This idea provides a short and straightforward proof. 



Lemma 1 (Equal angles) Let < <fi < ir. Let angle < d < tt/2 be so that sin($) = sin(0/2) sin(2#). Let 
m = — \\ and let Q' = Q(^4, Xi 0> 0)- Then 

\{y' x \Q,' m A\G)\ > l-a(2 + 4^ 2 m) 



(9) 



where = 

Let angle —it < p < ir be so that tan((,c/2) = tan(0/2)(l — 2a). Then \(f> — p\ < 2wa. Let Q = Q(-4, x, <t>i P) 
and let Tt denote the Hilbert space that Q and Q' act upon. Let || • || denote the operator norm on H defined by 
||0|| = sup{|0|r)| : | |r) | = 1} for any operator O on H. Then 



IQ'-Ql 





'1 




'1 















1 1 — e^-^l < 47r 2 a, 



(10) 



and thus ||Q' m - Q m || < 4vr 2 am. 

By definition of m we have that 1 1 - mi?| < f < 6, so sin(mtf) > VT~a and hence |(*i|Q m |*)| > 1 - 2a. Thus, 



K^'JQ'" 1 !*)! > K^IQ" 1 !*)! -47r 2 am > 1 - a(2 + 47r 2 m). 



(11) 



Lemma 1 follows. 

If we measure the state Q' m .4|0), then the outcome is good with probability at least 1 — 4a(l + 27r 2 m), which is 
at least 1 — 4a(^- + 11). The probability of measuring a bad state is thus upper bounded by 47r 3 ^ + 44a, which is 
(9(#) for < a < 1 and < 6 < tt. Theorem 2 follows. 
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Theorem 2 (Equal angles) Let A be any quantum algorithm that uses no measurements. Let a denote the success 
probability of A and let angle < 9 < -k/2 be so that sin 2 ((9) = a. Let tfi be any angle so that 9 < (f> < n. Let 
m = — |] where angle < i9 < n/2 is so that sin("d) = sin((^>/2) sin(2#). Let Q' = Q(.A, X, <j>, (/>)■ Then a 

measurement of Q' m A\0) will provide a good solution with probability 1 — 0(\t/o). 

Theorem 2 relies on two properties: Firstly, that operator Q(A, x> </>, (f>) approximates operator Q = Q(A, x, <f>, <f) 
sufficiently well. Secondly, that each application of Q implements a rotation by a sufficiently large angle d. Any 
set of angles (</>, ip) so that these two properties are fulfilled will provide a quantum amplitude amplification scheme. 
In general, the closer we pick the rotational angle to the maximal angle 29, the worse approximation for Q can 
we allow ourselves to use, and vice versa. In the next theorem, which is proven almost identically to Theorem 2, we 
express one way of capturing this duality. 

Theorem 3 (Any angles) Let A be any quantum algorithm that uses no measurements. Let a denote the success 
probability of A and let angle < 9 < tt/2 be so that sin 2 (#) = a. Let < <f> < tt and —it < if' < n be given angles. 

Let to = — 5] where angle < $ < n/2 is so that sm(i9) = sin(</>/2) sin(2#). Let — 7r < tp < tt be defined 
so that tan((p/2) = tan((/>/2)(l — 2a) and let 8 — \if' — if\. Let Q' = Q(A, x, <p, if') denote our approximation to 
Q(.A, Xi 0i <p)- Then a measurement of Q' -4|0) will provide a good solution with error probability at most 4a + e, 

provided <5<e 2 ^(^L) ^ 

The above theorem put bounds on the error 5 = \ip' — if\ that we can tolerate to still obtaining a quantum algorithm 
that succeeds with high probability. Suppose is a constant, say <f> = jn, then whenever 6 is at most Cy/a, for some 
appropriate constant c, the above algorithm finds a good solution with bounded error probability. Furthermore, if 6 
is at most 0(a), as it is if we pick if' — <p, then the error probability drops to being in 0(y/a). Thus, the smaller 
error in the choice of angles from the perfect case as expressed by Eq. (|B|), the smaller is the error probability of the 
algorithm. 

Theorem 3 is similar to the main result in [|| where it is proven that for quantum searching, the distance \(f> — ip\ 
must be at most of order for finding the marked element with constant success probability. Our result generalizes 
their result as it measures the error probability of the overall algorithm in terms of the distance from the perfect 
case in which Eq. (|) holds. The case 4> — f is already an approximation, which by itself introduces an error in 8 in 
the order of 0(a). In addition, Theorem 3 includes the cases that the angles <f> and ip are not constants but depend 
on a. However, the main benefit of Theorem 3 is that it is so easy to prove when one is given Theorem 1. Essentially, 
the proof of Theorem 3 reduces to bounding the distance between the two operators Q(^4, x> 4>i f) an d Q(^4, Xi 4>i f'): 
which is easy. 

VI. CONCLUSION 

Amplitude amplification is besides quantum Fourier transforms the most successfully used tool in quantum algo- 
rithms. It is a generalization of Grover's quantum searching algorithm and it allows a quadratic speed up of many 
algorithms. 

We prove that the phase condition in amplitude amplification can be expressed by the equation tan(i^/2) = 
tan(0/2)(l — 2a), where a denotes the success probability of the original algorithm. We show how to implement 
arbitrary rotations and how to boost quantum algorithms to succeed with certainty by utilizing angles (</>, if) that 
satisfy the above equation. In both cases, the number of iterations required increases linearly in the inverse of angle <f>. 
For instance, if we choose <j> = ~ for some constant c > 1, then we require ©(c-^) iterations to find a solution with 
certainty. 

Whenever the success probability a is not known a priori, we can approximate the above trigonometric equation 
by the linear equation (ft = if. This case has been studied by Long et al. |2|-|^] and in particular, they have shown that 
equal angles can be utilized in quantum searching. By considering the case cf> = if an approximation to the perfect 
case tan(yj/2) = tan(</>/2)(l — 2a), we can reprove the main results by Long et al. We believe that this approach 
makes the analysis and proofs straightforward and easy to understand. 
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